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■ Abstract. The Press-Schechter theory of the cosmological mass func- 
. tion and its modifications allow to constraint cosmological scenarios of 

the structure formation. Recently a few new models have been suggested 
that explored the influence of anisotropic collapse on the shape of the 
mass function. I discuss in more detail a particular model that assumes 
that a fluid particle becomes a part of a gravitationally bound halo when 
the smallest eigenvalue of the deformation tensor of the filtered initial 
density field reaches a certain threshold (like the filtered density contrast 
I reaches the threshold in the Press-Schechter formalism). Choosing the 

. smallest eigenvalue guarantees that the fluid particle in question experi- 

ences collapse along all three axes. The model shows a better agreement 

■ with the N-body simulations than the standard Press-Schechter model. 

o 
o 
o 

• 1. Introduction 

Q ■ The derivation of the distribution of masses in gravitationally bound objects is 

\ one of the princip le goals of the theory of the structure formation (for a review 

^ ' see onaco 1998| ). Comparing the theoretical mass function with observations 

provides important constraints on the cosmological models (see e.g. Bond &: 
Myers 1996, pahcall fc Fan 1998| , p.eichart et al. 19991 ). Rich clusters of galax- 



^ , ies represent a particular interesting class of objects for two reasons. Firstly, 

■ they are the largest gravitationally bound objects in the universe and therefore 

represent rare events. As one of the consequences of being rare events clusters 
are particularly sensitive to some parameters of the cosmological models 
and cTg). Secondly, the formation of clusters is relatively simple process since 
it is primarily determined by the gravitational dynamics while other processes 
(hydro, thermal, etc) are less important than e.g. in the process of galaxy for- 
mation. As a result the numerical simulations of clusters of galaxies are more 
realistic and reliable than simulations of galaxy formation. 

Measuring the mass function of galaxy clusters is not easy but recently a 
certain progress has been achieved for both optically (see e.g. Bahcall fc Cen| 



[1993| , Girardi et al. 199^ ) and X-ray ( pleiprich fc Bohringer 199£ ) selected sam- 



ples. Although there are systematic differences between mass functions obtained 
by different groups there is a general agreement in a broad sense. 

Most of theoretical derivations of the cosmological mass function are based 
on the ideas of Press & Schecter (1974) that can be summarized as follows: 
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The mass fraction ( F{> M) ) in gravitationally bound objects with masses 
greater than M can be estimated as the fraction of mass satisfied the 
collapse condition at this scale ( ^{6m > ^c) )■ F{> M) = 2^{6m > ^c) 



• The collapse condition is local i.e. it can be expressed in terms of the 
quantities at one point. 

• The quantity that determines the collapse is the linearly extrapolated fil- 
tered density contrast S]\f > 5c at a given point. 

• The value of the threshold 6c = 3/20(127r)^/'^ ~ 1.69 that corresponds 
to the collapse of a spherical top-hat model with the similar initial den- 
sity contrast. It was assumed that the collapse of the spherical top-hat 
model approximately corresponded to the virialization of the gravitation- 
ally bound clump. 

The mathematical aspects of the Press-Schechter formalism is outlined in the 
following section. Here I would like to discuss briefly some of ideas suggested 
since the formalism was proposed in 1974. 

The excursion set approach ( Peacock &: Heavens 1990 , Bond et al. 199i| ) 



justified the assumption that F{> M) = 2^{5m > <^c) in the case of a sharp 
A;-space filter. 

Many realized that the threshold 5c = 1.69 does not provide the best fit to 



N-body simulations. Although some authors used the canonical value (e.g. [Bond 
et al. 1991| , [Efstathiou et al. 1988|) others preferred the lower v alues: 5c = 1.58 
(Bond k, Myers 1996), 5c = 1.44 ^Carlberg fc Couchman 1989), or even as low 



as (5c = 1.33 ( [Efstathiou fc Rees 198^ and Klypin et al. 1995| ). Recently Shapiro 



et al. (1999) showed that the virialization of the top-hat model occurs when 
linear extrapolation of the density contrast reaches 5c ~ 1.52. 

One of the major efforts in reduction of the discrepancy of the theory with 
simulations has been related to incorporating the anisotropic character of grav- 
itational collapse. Bond & Myers (1996) developed a model that incorporated 



both the Zel'dovich approximation on large scales and the collapse of a homo- 
geneous ellipsoid on the nonlinear scale. Monaco (1995^ suggested a different 



collapse condition that corresponded to the collapse along the first axis in the 



Zel'dovich approximation. Audit et al. (1997) i ncorporated some of the non- 



linear effects into an anisotropic collapse model. Lee fc Shandarin (1998a) sug- 



gested to use the collapse condition corresponding to the collapse along all three 



axes as described by the extrapolation of the Zel'dovich approximation. Sheth 



fc Tormen (1999) obtained an analytic fit to the numerical mass function in the 
SCDM, OCDM and ACDM models and then gheth, Mo fc Tormen (1999] pro- 



vided a semianalytic derivation of the formula assuming an anisotropic collapse 
an incorporating some nonlocal effects. All but one models mentioned above 
assumed that the formation of a gravitationally bound object is related to the 



collapse along three axes. Only Monaco (1995) assumed the collapse condition 



corresponding to the collapse along only the first axis. 

In this talk I briefiy review the Press-Schechter formalism. Then I describe 
the derivation of the mass function (Aa-function) in the Zel'dovich approxima- 
tion. I compare th e result with the standard P ress-Schechter model and the 
model suggested by $heth, Mo k, Tormen (1999)] . I briefiy discuss the results of 
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comparison of the As-function with N-body simulations. Finally, I discuss the 
effect of the initial gravitational potential on the cosmological mass function and 
show that the clusters have a significant tendency to form in the troughs of the 
initial gravitational potential. 



2. The Press-Schechter Formalism 

The mass function n{M) is the number density of gravitationally bound clumps 
with masses between M and M + dM. Let F[> M) be the fraction of the 
mass contained in the gravitationally bound objects with masses greater than 
M. Press and Schechter (1974) suggested the fraction F(> M) and the mass 
function n{M) can be related as 

n(M) = — (1) 

where p is the mean mass density in the universe and the minus sign reflects the 
fact that is a decreasing function of M. 

Press and Schechter also made the assumption that the fraction of mass 
F{> M) can be estimated as a fraction of mass ^{6m > ^c) in the initial density 
field flltered with the window function W (corresponding to W in fc-space) 

6m{^, t) = D{t) J 6,ni^')W{\x' - x|/i?) d^x' (2) 

satisfying the collapse condition Sm > ^c- Here 6 = {p — p)/p is the density 
contrast, D(t) is linear growth factor, x is the comoving coordinate. The mass 
M and the linear scale R of the filter are related as 

M = fwY^^P^ (3) 

where fw is a factor depending on the shape of the smoothing fllter W. For a 
shar p k-space filter adopt ed here fw = 97r/2 w 14.1 and thus M = Gn'^R^p (see 
e.g. iLacey fc Cole (1994)| ). 



Assuming that the initial density contrast is Gaussian random field its pdf 
(probability distribution function) is 



p{^m) = ^= exp 

where the variance aj^j is a function of mass M 



(4) 



<yli = J j^P{k)W\kR), (5) 

where P{k) = \6k\'^ is the initial spectrum of perturbations and W{kR) is the 
window function in the k-space. 

Press and Schechter argued that a fluid element becomes a part of a gravita- 
tionally bound object of mass M when its linearly extrapolated density contrast 
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5m reaches the critical value 5c = 3/20(127r)^/^ ~ 1.69. This corresponds to the 
collapse of the top-hat spherical perturbation having the initial density contrast 
similar to the fluid element in question. The collapse of a spherical top-hat 
model has been assumed approximately to correspond the virialization of the 
clump. Recently Shapiro et al. (1999) showed that the virialization corresponds 
to (5c — 1-52 rather than to 5c — 1.69. 

The fraction of mass satisfying the collapse condition on scale M is 



^{8m > Sc 



27rcr(M) JSc 
5c 



exp 



2a2(M) 



d5M 



— erfc 
2 



^/2a{M) 



(6) 



where erfc(x) is the complementary error function. Assuming that F{M) 
"^{Sm > 6c) one easily obtains the mass function n(M) (eq. || ). 

One obvious problem with this result is that the normalization integral 



dF ^ ^{5 



M=oo 



>5c) 



(7) 



meaning that only a half of the mass is contained in the gravitationally bound 
clumps. Press and Schechter renormalized n{M) by introducing an additional 
factor of 2 (F(M) = 2^ {5m > 5c ) 



ripsiM) 



p p da 

MdM ~ '^MdM da 
2 p da 5c 



TiMdM a'^{M) 



exp 



51 



(8) 



Later, the normalization problem was correctly resolved in the frame of the 
excursion set model ( Peacock &: Heavens 1990| and Bond et al. 1991] ). The 
derivation of Press and Schechter did not take into account the so called cloud- 
in-cloud problem. Function ^'((5m > 5c) in eq.^ gives the fraction of mass that 
satisfies the collapse condition at the filtering scale M. However, some of the 
fiuid particles may satisfy the collapse condition at larger filtering scales. In 
the correct model the fiuid elements must be assigned to the clumps of mass 
Ml being equal to the largest filtering mass at which the collapse condition is 
fulfilled. In the excursion set formalism this corresponds to the first crossing of 
the collapse threshold 6c while 5 evolves with the growth of a. 

An elegant method to normalize the mass function was suggested by |Jedamzili 
(1995)1 ( see also the discussion in |Yano et al. (1996)D who derived the integral 



equation 



^{5m > 5c 



dMin{Mi)^P{M, Ml) 

M P 



(9) 



that relates the fraction of the fiuid elements firstly crossed the collapse threshold 
at the filtering scale Mi ( dMin{Mi)Mi / p ) and the fraction of mass satisfying 
the collapse condition at filtering scale M ( ^{5m > 5c) )■ Function P(M, Mi) 
is the probability that a fiuid particle firstly crossed the collapse threshold at 
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the scale Mi satisfies the collapse condition at the scale M. In the case of the 
sharp k-space filter and Gaussian 6m this probability is exactly equal to 1/2 for 
all Ml > M. Thus, the integral equation (^) can be immediately solved for the 
mass function n{M). The solution is the correctly normalized mass function of 
eq. I 



3. Mass Function in the Zel'dovich Approximation 

The simplest theory describing the anisotropic character of the gravitational 
collapse in a generic case of random initial condition is the Zel'dovich approxi- 
mation ( Zel'dovich 1970| , see also Shandarin &: Zel'dovich 1989 for a discussion). 



In particular, the Zel'dovich approximation provides a formula for an anisotropic 
collapse of a fluid element 

^ [1 - D{t)\M)][l - D(t)A2(q)][l - Z?(t)A3(q)] ' 

where D{t) is the linear growth function and Ai(q),A2(q) and A3(q) are the 
eigenvalues of the initial deformation tensor. Using the ordering convention 
Ai(q) > A2(q) and A2(q) > A3(q) the condition 1 — Z)(t)Ai(q) = has been 



interpreted as a collapse of a fluid particle along one principle axis ( Zel'dovich 
[L970| ). Similarly the conditions 1 — Z)(t)Aj(q) = (i = 2,3) can be interpreted 
as collapses along the second and third principle axes. 

yhandarm & Klypm (1984) showed that the formation of gravitationally 
bound clumps was the best correlated with the maxima of the smallest eigen- 
value (A3) of the initial deformation tensor. Although the formation of the 
clumps may be also related to other pointlike singularities ( Arnol'd et al. 1982| ) 



here we assume that a fluid particle becomes a part of a gravitationally bound 
clump of mass M when its smallest eigenvalue A3 reaches the critical value Ac at 
the largest filtering scale M (Lee & Shandarin 1998a). The Zel'dovich approxi- 
mation (eq.p!o|) predicts that the collapse condition is Ac = 1 (it is assumed that 
D{t) normalized to D{tQ) = 1, where is the present time). However, because 
of multistreaming effect all fluid particles (except the set of measure zero) enter 
the multi-stream flow regions before they collapse. We approximately incorpo- 
rate this complex effect by reducing the threshold Ac to a smaller value. The 
comparison with the Press-Schechter mass function as well as with the numerical 
mass function suggests that Ac = 0.37 is a reasonable choice. 

The derivation of the mass function in the Zel'dovich approximation is 
similar to the Press-Schechter derivation except that the collapse condition is 



A3(M) = Ac instead of 5m = 5c- Doroshkevich (1970) derived the joint pdf of 
three eigenvalues 

p(Ai, A2, A3) = exp - ^ + ^) (Al - A2)(A2 - A3)(Ai - A3), (11) 

where Ii = Ai -|- A2 -|- A3, I2 = A1A2 -|- A2A3 -|- A3A1 and cr^ is the density contrast 
variance as defined in eq. |5[ Integrating ^'(Ai, A2,A3) over two eigenvalues one 
can obtain the pdf of one of the eigenvalues. We are interested in the collapse 
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along the third axis and therefore ^(As) is of primary interest 



+ VS(20|-l)e.p(-g)erfc(V2^ 
- 20^Sexpf-M)|. (12) 



Repeating the derivation of the previous section using the pdf of eq.|T^ instead 
of eq.^ one arrives to an analog of the normalization integral equation (eq]9|) 

*(A3(M) > A,) = / dMin{Mi)^P{M,Mi), (13) 
Jm p 

here ^{\^{M) > Ac) is the fraction of mass where \^{M) > Ac on the filter scale 
M. 

Solving exactly eq.^ is much more difficult than eqj^ because now P{M, Ali) 
is not a constant. In the limit Mi — M <^ M the probability P = 0.5 as in eq.^ 
but in the hmit Mi > M it drops to P = 0.08. |Lee &: Shandarin (199"8a)| used 
the limiting value P = 0.08 and analytically derived the mass function in the 
Zel'dovich approximation 



n(M) 



25\/5 p da A 



3c 



247r MdMalj 



+ v/2^(2„|-l)»p(-||)e,.fc(v2^ 



A 



20^exp(-3)^. (14) 



a 



9A|\ 



2ay 



In the following sections I compare the obtained result with the Press- 
Schechter and Sheth-Mo-Tormen mass functions as well as with numerical sim- 
ulations. 



4. Comparison of Three Analytic Mass Functions 

Both the Press-Schechter and As-mass functions have a common factor -^-^ 
which depends on the initial spectrum and f{a) = dF/da that completely char- 
acterizes a model. Thus, comparing different models is convenient by comparing 
f{a) as functions of a. The Press-Schechter and As-functions are 



fps(.(^) = 




6 



25^/5 A 



3c 



247r cr2 



3 V Svr exp 



V 2a 



A 



+ \/2^ 20^-1 exp 



20— exp 




(15) 



Sheth and Tormen (1999) and Sheth, Mo and Tormen (1999) derived a new mass 
function that fits better the results of the N-body simulations 



/smt(o") = A 




2 ^/a6c 

exp 




(16) 



here a = 0.707, q 
condition 



0.3 and the constant A = 0.322 found from the normalization 



f{a)da = 1. 



(17) 



Choosing a = 1, g = and the constant A = 1 one obtains the Press-Schechter 
function. These three functions are shown in Fig. la. The small box shows the 
range of a where the theoretical mass functions were checked against N-body 
simulations by Sheth & Tormen (1999) and Lee &: Shandarin (1999)| . Fig. lb 
shows the ratios fps/fsMT and fx^/fsMT- 



5. Comparison with N-body Simulations 



Figure 2 shows the comparison of the Aa-function with the numerical mass func- 
tions for the scale invariant initial spectra: P{k) oc k"" with n = —1 and n = 
(see for the details Tormen (1998)). The top panel (n = — 1) shows a quite good 
agreement of the Aa-function with the numerical mass function, while in the 
n = case the agreement is significantly worse. Fig. 3 shows the comparison of 



the A3-fun ction with the N-body simulation of the SCDM model ( Governato et 
M- (1999)1) . At four epochs {z = 1.86,1.14,0.43, and 0) the As-function is in a 
better agreement that the Press-Schechter mass function. A similar result has 
been reported by pheth fc Tormen (1999) 



6. Large-Scale Biasing 

It has been noticed for sometime that the initial gravitational potential may no- 
ticeably affect the large scale structure. Kofman and Shandarin (1988) showed 
that the adhesion approximation predicts that the formation of voids is associ- 
ated with positive peaks of the primordial gravitational potential. Sahni et al. 
(1994) studied the effect and measured a significant correlation between the sizes 
of voids and the value of primordial gravitational potential in numerical simula- 
tions of the adhesion model. Recently, Madsen et al. (1998) have demonstrated 
by N-body simulations that the underdense and the overdense regions are closely 
linked to the regions with the positive and the negative gravitational potential 
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respectively. [Lee &: Shandarin (1998b) showed that the initial potential also 
affects the masses of clusters. 

In order to incorporate the primordial gravitational potential fluctuations 
term into the derivation of the mass function, we first derive the conditional 
probability density distribution p{S\ip < —ipc) {^c > 0): 



p{6\<p < -ipc) 



1 



2TTas 



exp 



1 + erf 




erf 



1 -1 



(18) 



Here aj, a^, 



and cr^ are the density, velocity and the potential variances re- 
spectively; K =< 6 ■ ip > /asfJip = (j'^/ascr^ is the crosscorrelation coefficient 
of the the density contrast 6 smoothed on the scale kc and the primordial (un- 
smoothed) potential fluctuations 99. As a result, eq|l| for the conditional mass 
function n{M\ip < —ipc) becomes 



OF das dF da^ \ 
d^sdM ^ doydM) 



(19) 



The further calculation needs to be done numerically. Fig. 4 illustrates how 
the mass function depends on the initial potential in the CDM model with 
r = = 0.25 normalized to erg = 1. The top panel shows the mass function 
for regions of positive and negative initial potential as well as unconditional 
mass function. The bottom panel show the ratio of conditional mass functions 
to unconditional one. 

We also calculate the probability that a clump with mass M is located in 
the potential regions satisfying the chosen condition, for instance, (p < —(pc 

pw < -^m = ^-iMk^^z^Pi^ < (20) 

where P{(p < —(fc) is the fraction of space satisfying the given condition (see 
Fig. 5). 

The scale of the initial potential 



\ 



Xdkk-^Pjk) 
' S^dkP{k) 



120/i~^Mpc 



f21) 



does not depend on any ad hoc scale; the dependence on ki is extremely weak 
(oc y/ln{l/ki) for the Harrison-Zel'dovich spectra assumed here). Th e geometry 



of the gravitational potential does not evolve much on large scales ( Kofman &: 
Shandarin (1988), [Madsen et al. (1998)| ). Therefore, the potential at present 



is very similar to the primordial one on scales much greater than the scale of 
nonlinearity. A simple explanation to this in the frame of the standard scenario 
of the structure formation is due to the fact that the mass has been displaced 
by the distance about 10/i~^Mpc ( ^handarin 1993| ). Therefore, the potential on 
scales greater than, say, 30/i~^Mpc has been changed very little. 
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For the model in question the scale of the primordial potential is found to 
be R^p ~ 120/i~^Mpc. The scale of the density contrast field reaches this value 
Rs = ^/^as/(Ts' ~ 120/i^^Mpc only after it is smoothed on kc ~ 0.017/iMpc~^. 
The corresponding density variance on this scale is iT5(0.017/iMpc^^) ~ 0.03. 
On the other hand, the number of clumps with masses 10^^ — lO^^/i~^M0 can 
easily be 30% greater in the troughs of the potential than the mean density 
n(> M) = 0.5[n(> M\ip < 0) + n(> M\ip > 0)] (see Fig. 4). Thus, the bias 
factor b (defined by the relation And/nd = bApm/pm) reaches at least 10 on 
the scale about 120/i~^Mpc. 

Figure 5 demonstrates that the most massive clusters (M > lO^^/i^^Mo) 
are almost certainly located in the the troughs in the initial potential. The 
bias defined as the density contrast of the clusters with respect to the mass 
density contrast b = 5ci/Sp reaches the value 3 — 10 on the scale of the potential 
w UOh-^Mpc (Lee & Shandarin 1998b). 



7. Summary 

In the talk I discussed new modifications of the Press-Schechter theory of the 
cosmological mass function. One assumes a different collapse condition that im- 
plies that a fluid particle becomes a part of a gravitationally bound object after it 
experiences collapses along three axes. The comparison with other models (Fig. 
1) shows that it predicts about 25% more gravitationally bound clumps than 
the Sheth-Mo-Tormen model in the range 0.45 > a > 3.1 where the comparison 
with the N-body simulations has been done. A direct comparison with the N- 
body simulations (Fig. 2 and 3) shows a quite good agreement although not as 
good as the Sheth-Mo-Tormen model. The Aa-function based on the Zel'dovich 
approximation has been obtained analytically similarly to the Press-Schechter 
function. A drawback of the derivation is a quite crude approximation of the 



probability function P{M,Mi) in the normalization integral eq.l3. A more ac- 
curate normalization will be reported separately. The Sheth-Mo-Tormen model 
also suffer from a normalization problem: the shape of the mass function has 
been derived but the normalization has been enforced by demanding equality of 
eq. 
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Another modification is the conditional mass function showing that the 
clusters of galaxies tend to form in the troughs of the initial gravitational po- 
tential and avoid the peaks of the potential. The gravitational potential field 
has a typical scale of about 120h~^ Mpc and as a result has an advantage of 
being independent of the arbitrariness of the smoothing scale (if the filter scale 
is smaller than roughly 50h~^ Mpc) and at present it has almost same geometry 
as at the epoch of decoupling. Figures 4 and 5 quantify this large-scale biasing. 
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Figure Captions 



Fig. 1. (a) The fraction of mass / = dF/da in the gravitationally bound 
objects as a function of a as predicted by the Press-Schechter model (short 
dashed line), Sheth-Mo-Tormen model (solid line), and As-model (long dashed 
line). The small box shows the range of a where the models have been checked 



against N-body simulations (see Fig. 2 in Sheth, Mo fc Tormen (1999) ) 
(b) The logarithm of the ratios fps/fsMT (short-dashed line) and fx-^/fsMT 
(long-dashed line). 

Fig. 2. The square dots represent the numerical mass function with poisso- 
nian error bars. The solid line is the Aa-mass function with Asc = 0.37 while the 
dashed, the dotted lines are the PS mass functions with 5c = 1.69, 1.5 respec- 
tively. The upper and the lower panels correspond to the n = —1 and the n = 



power- law models respectively. See also the top left panel of Fig. 2 in p?ormen 
(1998)1 . 



Fig. 3. The square dots represent the numerical data for the case of SCDM 
model with = 1, /i = 0.5. The solid line is our mass function with Asc = 0.37, 
and the dashed, the dotted lines are the PS mass functions with 5c = 1.69, 1.5 
respectively. 

Fig. 4. In the upper panel the conditional cumulative mass function sat- 
isfying chosen potential condition is plotted. The solid, the long dashed, the 
dot-dashed, and the dashed lines correspond to the conditions ip < —a^, (p < 0, 
if > 0, and if > respectively, while the dotted line represents the uncondi- 
tional cumulative PS mass function. The shaded area is la fit to the observa- 
tional cumulative mass function of rich clusters by Bahcall and Cen (1993). In 
the lower panel the ratio of the conditional cumulative mass functions to the un- 
conditional one is plotted for each condition. The CDM spectrum with F = 0.25 
normalized to as = 1 has been used. 



Fig. 5. The probability that a clump with mass M can be found in the 
regions satisfying chosen potential condition is plotted. The heavy solid, the 
heavy dashed, the solid, the dashed, the long dashed, and the dot-dashed lines 
correspond to the condition if < 0, if > 0, ip < —cr^, —a^ <ip<0,0<^p< a^, 
and (p > a^p respectively. 
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